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The aim of the studies to be made for this project is to determine 
the stress field in the lithosphere caused by the distribution of 
density anomalies associated to the geoidal undulations observed by the 
GEOS-3 and SEASAT Earth satellites in the Tonga region. 

Since the geoidal undulations do not determine uniquely the 
density-anomaly distribution causing the undulations-, different models 
of the lithosphere have been generated with different assumptions on the 
density distribution and geometry, all generating a geoid profile almost 
identical to the observed one. 

The first model used is that with the Airy isostatlc hypothesis 
(Christian 1984), it consists of a crust (defined as a lower-density 
layer) of density 2.85 laying on a lithosphere of density 3.35. The 
models obtained with different compensation depths give residual 
shortwavelength anomalies of the order of several tens of mgal and 
several tens of meters geoidal undulations. This clearly indicates that 
in the Tonga region there is no isostasy of the Airy type because the 
observed geoid has very smooth undulation of about 25 m over a distance 
of 2000 km. 

We also used the Pratt isostatic hypothesis in a model consisting 
of a crust of variable density laying on a lithosphere of higher 
density. This model gives smaller residual anomalies (Wainright 1983) 
but still shows that in the Tonga region there is no isostasy of the 
Pratt type because the observed geoidal undulation are much smaller and 
smoother than the residual undulations associated to the Pratt model of 
isostasy. 
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It is thus clear that the density anomalies needed to produce the 
observed geoid are not producing isostatic equilibrium at any depth, and 
therefore there is no hydrostatic equilibrium. The estimate of the 
consequent deviatoric stress is the aim of the research. 

To compute it one must therefore compute models of the lithosphere 
which reproduce the geoidal undulations observed by GEOS-3 and SEASAT. 

The first model (a) (Fig.l) was generated by assuming that the 
density of the crust is 2.85 (Carlson and Raskin 1984) and that of the 
lithosphere below it is 3.35 and taking into account the density (or 
thickness) variation with time of the oceanic lithosphere (McAdoo and 
Martin 1984) . This model leads to an oceanic crust 24 km thick landward 
of the trench and 32 km thick landward of the trench. 

In a second model (b) (Fig. 2) the density of the crust is assumed 
2.85 and that of the lithosphere below 3.50; this model gives an average 
thickness of 18 km for the crust seaward of the trench and 28 km 
landward of the trench. 

This raises serious questions on the meaning of the discontinuity 
surface in the velocities of the seismic wave found at the depth of 6 km 
below the ocean bottom (Shor et al. 1970, Turcotte and Shubert 1982) and 
assumed to define the thickness of the crust; it is clear now that a 
better definition of the term crust is needed in order to distinguish 
between the discontinuities of the different physical parameters. 

In all cases examined the oceanic lithosphere (Fig. 3, 4) has been 

“3 3 

assumed to vary in density with time at a rata of. about 0.2 10 gr/cm 
Myears (rather than in thickness) and this low rate density variation 
together with the thinner crust landward of the trench and the thinner 
downgoing slab are responsible for the long wave length geoidal 
undulation observed by GEOS-3 and SEASAT. 
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In both cases (a) and (b) the computed models of the geoid deviates 
from the observed by less than 2 meters . 

The density distribution of models (a) and (b) have then been used 
to determine the maximtim shear stress (mms) field caused by their load 
on the layers below. 

To estimate the stress field in a layered spherical Earth model we 
used the method of Caputo (1961); however the formulae of that paper 
would not allow for the body force associated to buoyant masses and had 
to be modified accordingly (Caputo 1984) ; the new formulae allow to 
compute the stress field in a layered sphere caused by the most general 
distribution of surface tractions and body forces. 

The most recent results (Caputo et al. 1984) obtained using these 
formulae indicate that the mss caused by a mountain range and by its 
isostatic compensation, when this is exactly one below the range, is 
limited to the layer under the load including the layer containing the 
isostatic adjustment (Tables 1, 2); when the isostatic mass is displaced 
with respect to the load then the mss extends to a much greater depth 
(Tables 3, 4). In all cases the mss is at the most one third of the 
load. It is to be noted that when the isostatic mass is exactly below 
the load the mss reaches its maximum at the surface of separation 
between the top layer and that containing the isostatlc mass and then it 
decreases linearly to almost zero within the latter layer (Fig. 6). 

We then applied these results to the density anomaly models across 
the Tonga Trench obtaining a first estimate of the mss at several depths 
below the sea surface. To do this we computed, at depth of 9 km, 31 km, 
59, 80 km, the pressure caused by the masses above the respective depth 
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for all the density models previously obtained for the Tonga region. 

In all cases we found that the most important features (Fig. 5) were 
pressure variations over relatively short horizontal distances which 
could be simulated by a box-like normal-traction (positive or negative) 
as considered by Caputo et al (1984) or a step like variation. The 
latter being the less relevant case in terms of the amount of traction 
variation. 

The mss found in the Tonga region has several features. Maximum 
shear stress (mss) of about 300 bar are found close to the Trench 
(region A) , and in the back area region (region B) extending to a great 
depth in the volume between them. 

The analysis of the seismic activity in the Tonga region reveals 
that this is concentrated in regions A and B. In regions A and B there 
is normal depth seismic activity with moderate to large magnitudes , 
while between the two regions, as well known, the seismicity extends to 
a depth about 650 km. 

It is important to note that there are other places in the Tonga 
region where the mss reaches values of 200 bar and there is no evidence 
of occurrence of large earthquakes in recent times namely in the area 
about 800 km landward of the trench. Detailed analysis of the 
seismicity of these regions may possibly reveal that they can be 
considered as locations of possible seismic gaps. 

In this research it has also been found that the portion of the 
lithosphere under the seamount, in the cross section analysed, is not in 
isostatlc equilibrium, and that the mss under the seamount is of the 
order of 250 bar at a depth below 50 km. 
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The stress field under seamounts had been estimated by Lambeck and 
Nakibogln (1980) and previously by other authors who modeled the 
deflection of the ocean lithosphere and gave values of more than 10 kbar 
for the mss (Walcott 1970, 1976, Watts et al. 1975); however Lambeck and 
Nakibogln (1980) argued that there are some plausible methods of 
reducing this mss and suggested that the mss should not exceed 1 kbar. 
The hypotheses invoked by Lambeck are lower density for the sediments 
fill in, depth-dependent nonelastic rheology and large-deflection, 
theory for the large loads. 

The models determined in this research confirm that even without 
invoking the causes hypothesized by Lambeck and Nakibogln (i960) the mss 
under the seamount should always be about one third the load caused by 
the mount; which in the case of the seamounts studied by Lambeck would 
imply a mss less than 0,5 kbar. 

It well known that rheology may set constraints on the mss that may 
have been accumulated in a region. It also known that there is no 
evidence that the Mantle and/or the Lithosphere may be described as 
Maxwell or Standard Linear Solids although several authors have 
tentatively used these models in their studies. In order to determine 
the effect of’ other possible rheologies we studied models different from 
the Maxwell and the Standard Linear Solid (Caputo 1984a); some of these 
models are supported by laboratory data (Caputo 1984b, Caputo 1984c). 

Concerning the effect of rheology in the estimate of the mss of the 
lithosphere we must note that the theoretical research conducted in this 
project (Caputo 1984a) has shown that the relaxation time t: defined as 
the time to reduce the stress to e ^ of its initial value is not:. 


5 


indicative of the state of stress of a material at times t >> t > 
because there are rheological models which would maintain the e ^ 
reduced value of the initial stress for an almost Indefinite time (Fig. 
7> Caputo 1984a). Therefore the stress estimates computed in this 
research may have to be increased due to the residual stresses 
accumulated in previous time. It is also implied (Caputo, 1984d) that 
other phenomena such as the elastic rebound and their effect on the 
may have time history different from that assumed; for instance the 
variation of the J 2 may be due to the superposition of the effects of 
more than one glacial period. 

The research made for this project concerning the rheological 
models of the mantle has also show that using the rheological models 
resulting from laboratory data on granite the response of the mantle to 
applied stress fields is such that relascation times is almost 
Independent of the wavelength of the stress field applied (Caputo, 
1984b, Fig. 8). 

A final note on the subductlon zone process is due. Giardini and 
Woodhouse (1984) discussed 17 moment tensor solutions in the Tonga 
region and their implications for deformation within the subductlon 
zone; they found a complex cross-cutting pattern of interacting shear 
bands and concluded that the Benloff zone should be seen "as that part 
of the convective flow which by virtue of temperature, composition and 
strain rate accomplishes its deformation through episodes of shear 
instability. " 

After the findings of this research, namely the theoretical 
estimates of the mss under mountain ranges (Caputo et al. 1984), the 
models of the lithosphere in the Tonga region and the mss distribution 
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underneath we may comment that the idea of the slab smoothly bending 
with very slowly varying density and or thickness defined by the 
seismicity distributions may need a revision but in a direction 
different than that of Giardini and Woodhouse (1984), One could view 
part of the seismicity as caused by the density anomalies distributions 
through the mss generated in the lithosphere and in the deeper layers, 
as the distribution of the latter is largely superimposed over the 
volume where the former occur. 

Better evidence on this hypothesis could be supplemented by a 
detailed modelling of the mss under the Tonga region. This must now be 
obtained by computing it directly for the models of the density 
distribution of the lithosphere already determined rather than using the 
results of computations made for schematic models. More evidence should 
also be found by testing this hypothesis in other subductlon regions of 
the world. 

At this stage we have already obtained the density anomalies in two 
more profiles perpendicular to the Tonga Trench and we are in the 
process of computing the load at various depths due to these anomalies. 
From the load we shall infer the mss at depth. 

The computer programme for the calculation of the mss in a layered 
spherical Earth has been completed to Include the body forces in the two 
upper layers. This will allow to study more complex problem an obtain 
more realistic results. 
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The Research associated to the solution of the problems of this 
project led to the publication of the following papers; 

Topography and its isostatic ccqpensation as a cause of 

seismicii^. Tectonophysics , 79.. 73-83, 1983, Caputo M, , 
Milana G, , Rayhom J, 

Relaxation and free modes of a self-gravitating planet. 

Geophys. J. R. ast. Soc. , 77, 789-808, 1984a. M. Caputo. 
Topography and its isostatic corpensation as a cause of 

seismicity, a revision, Tectonophysics, 111, 25-39, 1985a. 
Caputo M. , Manzetti V. , Nioelli R. 

Spectral Rheology, Proceeding Synposiim Space Techniques for 
Geodynamics, Hungarian Acad, of Sc., J. Scmogi, 

C. Rigberg Ed. , 1984b. M. Caputo. 

Determination of the creep, fatigue and activation energy 

fron constant strain rate experiment. Tectonophysics, 91, 
157-164, 1983, M. Caputo. 

Nonlinear and inverse rheology of rocks, 1985c. M. Caputo, (In 
press) . 

Generalized rheology and geophysical consequences, Tectonophysics, 
1985d. 

and to the following theses: 

Wainright E.J. , Calculation of isostatic gravity anomalies and geoid 
heights using two dimensional filtering; inplications for 
structure in subduction zones, Ph.D. Dissertation, Graduate .College 
of Texas A&M University, 1983. 
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Christian B.E., Utilization of satellite geoldal anomalies for computer 
analysis of Geophysical mode^j. Master thesis. Graduate College of 
Texas A&M University, 1985. 

Sevuklekln M.T., The estimate of the shear stress field of a self 

gravitating elastic planet during its formation. Master Thesis, 
Graduate College of Texas A&M University, 1984. 

Martin Robert, The propagation of the stress due to the ridge push. (In 
progress) . 

Mecham Brent, The distribution of maximum shear stress is subduction 

zones due to topography and isostatic compensation. (In progress) . 
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OR POOR QUALl't^V 


TABLE 1 

Maxim uaj shear stress to the spherically layered Earth model described in (he (ext. The fint and the 
second layers are 30 km thick. The load of 330 bar is applied (o a strip of the outer surface defined by lac 
buoyant mass in the second layer Is actins between the same latitudes. The contoun 
are at 20 bar, 70 bar, 120 bar. 
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OF POOR QUALITY’’ 


TABLE 3 

Muiattffl ibcv nrcu in (be sphericailJy Uyocd Eutb model described in the tezL The fim and second 
Uytn are M km tiddc, Tbe lo^ of 3S0 bw is applied to a strip of tbs outer surface defined by the Jat 
^is^ind Tbe buoyant mass in (be teeoad layer is actins between the latitudes and 


Tbe oosiouts are at 20 bar, 70 bar. 120 bar. 
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Top layer: 

X - 3.8 • 10”. ;i - 3.8 • 10” , r , « 6.371 • 10*, - 6.341 • 10* 

Second layer: 

X - 7.5 • 10”, 11 - 6 A- 10”, r, - tutriable, - 6.341 ■ 10* 
Core; 

X “ 6 ■ 10*^, p.^6’ 10*^, T] » variable, rg ■= 0 


TABLE 4 


As in Table? nith the second layer 80 km thick, syminein'c case. 













waBrnHinSOr*— 3^ 4flr*>D—r“r'"'35 





co3iin*iinSO0--?q oi->C>— ri — S cnsmHmS 









coajmHms 





Figure 4 • Ccanplate Slab-: 
residual geoid. 
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Figure S- . Pressure at depths o£ 9, 31, 55, and 80 kilometers. 






OF POOR QUALllY 



Fig 6 Alanxsuin ibear nnu u huciioa of depth in the seoand layer of the model (beginmng at 30 bn 
depiny gtven in Fig. 1. In eaoh of tho four cases considered for the thiclcness of the second layer, the depth 
measured front ihe'top of the second layer is normalized to the thickness of ibe layer. Full circles refer to 
a second layer with 30 km thickness j squares refer to 55 km thjdcaess j triangles refer 

to SO km ihicbicss ^ and open circles refer to 175 km thickness Ai Tru bat^» 
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Figure 7. Reiaxafcion time f(t) of the rheological models identified by the 

parameter Z; the relaxation time is assumed unit and is arbitrary; here 
ia the same for all Z, only the time history of the relaxation 
cepends on Z. 



